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ABSTRACT 

We present a detailed investigation into the recent star formation histories of 5,697 
Luminous Red Galaxies (LRGs) based on the H<5(410lA), and [O n] (3727A) lines and 
the D4000 index. LRGs are luminous (L> 3L*), galaxies which have been selected to 
have photometric properties consistent with an old, passively evolving stellar popula- 
tion. For this study we utilise LRGs from the recently completed 2dF-SDSS LRG and 
QSO survey (2SLAQ). Equivalent widths of the H<5 and [O n] lines are measured and 
used to define three spectral types, those with only strong EW absorption (k+a), those 
with strong [O n] in emission (em) and those with both (em+a). All other LRGs are 
considered to have passive star formation histories. The vast majority of LRGs are 
found to be passive (^80 per cent), however significant numbers of k+a (2.7 per cent), 
em+a (1.2 per cent) and em LRGs (8.6 per cent) are identified. An investigation into 
the redshift dependence of the fractions is also performed. A sample of SDSS MAIN 
galaxies with colours and luminosities consistent with the 2SLAQ LRGs is selected to 
provide a low redshift comparison. While the em and em+a fractions are consistent 
with the low redshift SDSS sample, the fraction of k+a LRGs is found to increase 
significantly with redshift. This result is interpreted as an indication of an increasing 
amount of recent star formation activity in LRGs with redshift. By considering the 
expected life time of the k+a phase, the number of LRGs which will undergo a k+a 
phase can be estimated. A crude comparison of this estimate with the predictions 
from semi-analytic models of galaxy formation shows that the predicted level of k+a 
and em+a activity is not sufficient to reconcile the predicted mass growth for massive 
early- types in a hierarchical merging scenario. 
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1 INTRODUCTION 

The widely accepted paradigm for galaxy formation is the 
hierarchical merger model which dictates that galaxies form 
via continuous mergers with smaller objects (e.g. Cole et al. 
1994). This scenario is a natural consequence of the A cold 
dark matter (ACDM) model of structure formation, which 
has been highly successful in explaining both the observed 
structure from the large redshift surveys (2dF Galaxy Red- 
shift Survey (2dFGRS); Colless et al. 2001, Sloan Digital Sky 
Survey (SDSS); Stoughton et al. 2002), and the anisotropy 
in the cosmic microwave background (Spergel et al. 2003). In 
the context of this model one would naively expect the most 
massive galaxies to also be the youngest as they have been 
'built up' slowly via mergers. However traditionally early- 
types are thought to have formed early and rapidly, with 
a vast amount of observational evidence points supporting 
this conclusion. The homogeneous nature of massive early- 
type galaxy properties, such as small scatter in the colour- 
magnitude relation (CMR) (Visvanathan & Sandage 1977; 
Bower, Lucey & Ellis 1992; Terlevich, Caldwell, & Bower 
2001; Aragon-Salamanca et al. 1993; Stanford, Eisenhardt, 
& Dickinson 1995; Ellis et al. 1997; Stanford, Eisenhardt, & 
Dickinson 1998, among others) and fundamental plane (FP) 
(Dressier et al. 1987; Burstein et al. 1997), and the lack 
of evolution in these relations with redshift (Kodama et al. 
1998; van Dokkum & Stanford 2003), is consistent with a 
scenario in which most, if not all massive early-types, formed 
at high redshift (at least z > 2) and have experienced only 
passive evolution of their stellar populations since. 

While this could simply be a result of selection bias 
(van Dokkum & Franx 2001), the lack of direct evidence 
of strong evolution in the population of early-type galax- 
ies since z ~1 remains a serious problem for the hierarchi- 
cal model of galaxy formation. In recent years much work 
has gone into reconciling these differences. While theoretical 
models have moved towards consistency with the CMR and 
FP (Kauffmann & Chariot 1998; De Lucia et al. 2005, among 
others), observational studies based on the spectroscopic and 
morphological properties of early-types have found signifi- 
cant populations which show evidence for recent star forma- 
tion and/or interactions. A morphologically based study by 
Michard & Prugniel (2004) found a significant number (~ 30 
per cent) of nearby (z ~ 0.01) early-type galaxies showed ev- 
idence of recent interactions/mergers, although 40 per cent 
of these showed no evidence for a young stellar population. 
A similar study of 86 nearby bulge-dominated red galaxies 
by van Dokkum et al. (2005) found 71 per cent had evi- 
dence of tidal interactions. However given that both of these 
studies were undertaken at low redshift the volume sampled 
was too small for significant numbers of LRG analogues to 
be present. Moving to higher redshifts, Willis et al. (2002) 
found evidence for current star formation, in the form of de- 
tection of the [O n] emission line, in over 20 per cent of a 
sample of 415 luminous field early-types at z ~ 0.3. While 
these studies show that evolution is occurring in at least 
some of the early-type population, it remains unclear if all 
early- types have experienced evolution of this type. Another 
contentious point is whether early-type galaxy acquire mass 
predominately in the form of star forming gas or already 
formed stars via so called 'dry mergers' ( Bell et al. 2005). 

In this study we look for evidence of recent star forma- 



tion in a sample of Luminous Red Galaxies (LRGs). LRGs 
are luminous (L > 3L*) early-type galaxies, analogous to 
bright cluster galaxies (BCG), selected via their red colours 
which are broadly consistent with a passively evolving stel- 
lar population. These galaxies are ideal candidates for the 
purposes of this study as they have been shown to have rel- 
atively homogeneous spectral properties (Eisenstein et al. 
2001), enabling any deviation due to recent star formation 
episodes to be readily identified. They are also the most mas- 
sive galaxies, which, if the hierarchical merger model is to be 
correct, must also have the most vigorous merger histories. 
Thus a sufficiently large sample of LRGs provides the perfect 
test bed for current theories of galaxy formation/evolution. 

Here we probe the star formation histories of our 
LRG sample via their spectral properties, in particular the 
[O ll](3727A) emission and H<5(410lA) absorption lines and 
the D4000 index. LRGs are selected from both the SDSS 
and 2dF-SDSS LRG and QSO (2SLAQ) surveys, result- 
ing in a sample of 5,697 LRGs covering a redshift range of 
0.2 < z < 0.7. Details of the sample selection are presented 
in Section |21 Section |3] describes the analysis and modeling 
performed, with the subsequent results presented in Section 
2] The implications and conclusions to be drawn from these 
results are then presented in Sections |S] and|5] respectively. 

Where necessary, we have assumed a flat cosmology 
with n m = 0.3, Q A = 0.7 and H = 70 kms^Mpc" 1 . 



2 THE 2SLAQ LRG SURVEY 

Completed in August 2005, the 2SLAQ survey has col- 
lected spectra and measured redshifts for over 10,000 high 
redshift QSOs and ~ 14,000 LRGs in the redshift range 
0.45 < z < 0.8. The survey is essentially a high redshift 
extension to the SDSS LRG survey, utilising the 2dF fa- 
cility on the Anglo-Australian Telescope to go deeper than 
was possible with the SDSS survey telescope. The 2SLAQ 
dataset has already been utilised to produce photometric 
redshifts (Padmanabhan et al. 2005), and a large photo-z 
LRG catalogue (Lahav et al. 2006), as well as probe the 
evolution of the LRG luminosity function from z = 0.2 to 
z = 0.55 (Wake et al. 2006). 

While a detailed description of the survey is presented in 
Cannon et al. (2006), we briefly review the pertinent points 
here. The 2SLAQ survey was designed to be an extension 
of the SDSS LRG survey which sampled LRGs in the red- 
shift range 0.15 < z < 0.5, to higher redshifts (Eisenstein 
et al. 2001). The SDSS LRG survey utilised two colour se- 
lections in g — r and r — i colour space. One colour selection 
was used to sample LRGS at redshifts less than 0.4 (cut I), 
where the 4000 A break features prominently in the g band, 
while another (cut II) was used for higher redshifts, where 
the 4000A break is located in the r band. The two cuts were 
necessary as the 4000A break dominates the colour evolu- 
tion of massive early-type galaxies, such as LRGs, and as 
such the transition from one band to another causes the 
track taken by the galaxy in colour-colour space to change 
considerably. In addition to the colour cut, a sliding mag- 
nitude limit was utilised for the low redshift cut I LRGs as 
the colours of LRGs at z < 0.4 are quite similar to less lu- 
minous star forming galaxies at lower redshifts. The 2SLAQ 
survey utilised a colour selection similar to the higher red- 
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shift, cut II, SDSS colour cut with some key differences. As 
the number of fibres available to 2SLAQ survey per sq deg. 
was significantly different (~ 67 per sq deg. for 2SLAQ as 
compared to only ~ 2 per sq deg. for SDSS cut II targets) 
the colour selection has to be 'loosened' to provide enough 
targets to fill a 2dF field. This can be seen in FigureQwhich 
shows the position of the 2SLAQ LRGs used in this study 
with respect to both the 2SLAQ and SDSS cut II colour cuts. 
A passive evolution track, provided by Bruzual and Chariot 
(2003) models, is shown as the solid dashed line. The SDSS 
LRG survey, with its smaller number of targets per field cuts 
very close to the passive evolution track and thus goes close 
to only selecting 'red and dead', passively evolving, galaxies. 
By contrast the increased availability of fibres from 2dF al- 
lows the 2SLAQ selection to be quite broad in g — r colour. 
This has a distinct advantage over SDSS in that LRGs with 
a range of star formation histories are allowed by the 2SLAQ 
selection, not just strictly passive evolution. In addition to 
the colour cuts a magnitude limit of ioev < 19.8 is imposed 
so as to ensure that reasonable signal to noise (~ 3) could 
be achieved in 4 hrs. 

Another key difference between the 2SLAQ and SDSS 
cut II selections is the grouping and prioritization of targets 
based on a number of colour cuts. LRGs which lie above the 
line marked 'Sample 8' were the main goal of the 2SLAQ 
survey and given the highest priority in targeting. The jump 
in r — i from the SDSS cut II selection and the 2SLAQ 
sample 8 cut is a result of the desire to acquire distinctly 
higher redshift LRGs than those in SDSS. Sample 9 LRGs 
were given lower priority and are essentially lower redshift 
targets closer to the redshift range of the SDSS cut II LRGs. 
As sufficient sample 8 and 9 targets were not necessarily 
available in each 2SLAQ field a small number of 'fibre-fillers' 
were selected below the Sample 9 cut. These are even lower 
redshift LRGs in the region sampled by SDSS cut I. 

The 2SLAQ observations were performed using exclu- 
sively the 2dF spectrograph 2 with a 600 lines per mm V 
grating. This gives a dispersion of 2.2A per pixel and a res- 
olution of about 5A. The detector used was a Tekl024 CCD 
with 1024 x 1024 pixels. Thus using 6150A central wave- 
length the resulting spectra have a wavelength coverage from 
5050A to 7250 A. This coverage was chosen to ensure that 
the Ca H&K lines were always present in the spectra across 
the target redshift range of the survey. 

In this study we select 2SLAQ LRGs which fall into the 
Sample 8 and Sample 9 selections. In addition only LRGs 
which have confident redshifts are selected by requiring the 
redshift quality flag be greater than or equal to 3. Quality 
flags on 2SLAQ LRGs are set by the redshift code Zcode, 
a derivative of the code used to determine redshifts for the 
2dFGRS survey. A quality flag of 3 or greater indicates a 
high level of confidence in the stated redshift (see Cannon 
et al. 2006 for details). 

We also require that both the [O n] and H5 lines are 
within the wavelength coverage of the spectra and that nei- 
ther line is in close proximity to the strong sky emission lines 
at 5577A, 5900A, 6300A and 7244A. The signal-to-noise ra- 
tio per pixel, averaged across the spectrum, is also required 
to exceed 3.5 as it was found that lower signal to noise data 
produces too many false line detections, as will be discussed 
in Section |4~T1 

These requirements bring the number of 2SLAQ LRGs 
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Figure 1. Colour-colour diagram for LRGs selected in this study. 
The solid lines show the region selected by the SDSS cut II, the 
dashed lines represent the 2SLAQ selections. The dot-dashed line 
shows the track taken by a passively evolving early-type galaxy 
model. The circles are located at Az = 0.1 intervals, starting at 
2 = 0.2 (lower left) and ending at 2 = 0.6 (top right). The sharp 
bend in the track near 2 = 0.4 occurs where the 4000A break goes 
from the g to the r band. 

utilised in this study to 5,697. It should be noted that 
all magnitudes and colours throughout are on the ABmag 
scale, in line with the 5-band SDSS photometric system (see 
Fukugita et al. 1996). 



3 METHODOLOGY 

3.1 Measuring Equivalent Widths 

Equivalent widths (EW) of known spectral features are mea- 
sured using a pseudo-Lick passband flux adding approach. 
The flux in a set band around the known line wavelength is 
computed and compared to a continuum level. The contin- 
uum is derived by calculating the average flux level in two 
sidebands and interpolating across the region of interest us- 
ing these values. For the H<5 and H7 lines we use the Lick 
standard H<5a and H7A definition for the passband and side- 
bands (Worthey & Ottaviani 1997), while for the [O 11] line 
we use the definition from Balogh et al. (1999). These defi- 
nitions are presented in Tabled It should be noted that the 
passband definitions for the absorption lines are the same as 
those used in the output from Bruzual and Chariot (2003) 
models. Note that throughout this paper we define a neg- 
ative equivalent width as emission and a positive value as 
absorption in a stated line. 

D4000 break strengths are also measured for each spec- 
trum in the sample. We define the D4000 index as per Balogh 
et al. (1999), namely the ratio of the flux in two 100A win- 
dows centred on 4050 and 3900 A. While the 2SLAQ spectra 
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are not flux calibrated, we find that the change in spec- 
trograph response is insignificant compared to the intrinsic 
error in the spectrum and thus the D4000 indices are unaf- 
fected. 

Errors on the equivalent widths are calculated utilis- 
ing equation A8 from Bohlin et al. (1983). Errors in the 
continuum are based on the standard error in the mean of 
the sideband flux; errors in the line are based on variance 
spectrum calculated by the 2dFDR pipeline. It should be 
noted that while the equivalent width, and D4000, errors 
take into account 'random' errors such as photon counting 
statistics, they do not include contributions from systematic 
errors in the spectra, such as poor sky-subtraction or flat- 
fielding. Contributions from these sources to the equivalent 
width measure could be as high as a few percent, and as 
such the quoted errors are likely underestimates of the true 
error. While this may be of concern for future users of the 
data presented in Table 1, this is not considered a problem 
here as the error on the indices does not play a major role 
in any of the results presented. 

A listing of the equivalent width measures on the 
[O n] and H<5 line, as well as the D4000 indices for 
all 5,697 2SLAQ LRGs used in this study can be at 
http : / /lrg. physics .uq. edu. au/publi .html 



4 RESULTS 

4.1 Definitions based on [O iljand H<5 

In this section we present a number of spectral classifica- 
tions based on the [O n] and H<5 equivalent widths which 
will be used in the subsequent analysis. We split the LRGs 
into 4 categories based on their position in the EW{H<5}- 
EW{[0 n]} space. 

(i) Passive - LRGs with no significant [O n] emission 
(EW{[0 ll]}> -8A) or US absorption (EW{H5}< 2A), in- 
dicating an old passively evolving stellar population. 

(ii) k+a - LRGs with significant H<5 absorption 
(EW{FM}> 2A) but no [O n] emission (EW{[0 n]}> -8A). 

(iii) em - LRGs with significant [O n] emission 
(EW{[0 n] }< -8A) but no U8 absorption (EW{H5}< 2A). 

(iv) em+a - LRGs with significant H<5 absorp- 
tion (EW{H<5}> 2A) and significant [O n] emission 
(EW{[0 n]}> 2A). 

The implementation of this classification scheme quan- 
titatively involves a careful determinination of what a 'sig- 
nificant' amount of [O n] emission or H<5 absorption is in 
light of the quality of the dataset. Physically the divisions 
in EW{H<5}-EW{[0 n]} space should be designed in such a 
way to isolate LRGs which have star formation histories that 
are distinctly different from simply passive evolution. In the 
high signal-to-noise, and resolution limit this is relatively 
easy as passive evolution models have no [O n] emission and 
low levels of H5. However in the low signal-to-noise, low res- 
olution regime, in which the 2SLAQ LRG spectra predom- 
inately lie, variations in the EW{H5} and EW{[0 n]} will 
be dominated by noise rather than intrinsic physical prop- 
erties. Thus rather than use a physical justification for our 
EW{H5} and EW{[0 n]} divisions, we calculate threshold 
values at which it is unlikely that an LRG has a spectrum 



of a passively evolving stellar population in the presence of 
noise. 

To achieve this we introduce 'mock' H<5 and [O n] pass- 
bands into our equivalent width analysis. The wavelengths 
of these 'mock' passbands are positioned in quiet regions of 
the LRG spectrum away from known spectral lines such that 
any significant equivalent width measures will be purely a 
result of noise. The passbands and sidebands used in the 
mock analysis are found in Table Q For each spectrum in 
the 2SLAQ LRG sample the equivalent width of the 'mock' 
line index is measured using the method described in Sec- 
tion 13.11 As well as giving us a benchmark with which to 
determine the optimum EW{Hi5} and EW{[0 n]} threshold 
values, this analysis has the advantage that for any arbi- 
trary set of EW{H5} and EW{[0 n]} thresholds we have a 
good estimate for the false detection rate, i.e. the number of 
LRGs we expect to fall into our k+a, em+a or em sample 
due simply to noise. 

Figure [5] shows the distribution of [O n] and H<5 equiv- 
alent widths for both the real and mock index definitions. 
Note again that negative equivalent widths indicate emis- 
sion while positive indicate absorption. The low signal to 
noise of the 2SLAQ data is apparent in the spread of the 
mock equivalent widths for both lines. The distribution of 
mock equivalent widths covers the same range as the real 
equivalent widths for both H<5 and [O n], meaning it is im- 
possible to define any divisions which will eliminate all false 
positives. Thus we try and define divisions that will min- 
imise the number of false positives while still selecting a 
statistically significant number of LRGs in each of our clas- 
sifications. It can be seen that the H<5 mock equivalent width 
distribution drops significantly at ~ 2A, creating a large dif- 
ference between the real and mock distributions. If we use 
this value to define 'significant' H<5 absorption the number 
of combined false positives in the k+a and em+a categories 
is 102, while the number of LRGs with real EW{H5}> 2A 
is 330. Thus using a division at EW{H<5}=2A results in 30.1 
per cent of k+a and em+a's being false positives. If we de- 
crease the division to lA the ratio of false positives to real 
detections increases to 33.5 per cent, while increasing the 
threshold to 5A only decreases the ratio to 29.9 per cent, 
hardly worth a decrease in number of real detections of 275 
(600 per cent). 

By a similar justification we define the EW{[0 n]} di- 
vision. The distribution of mock Oil equivalent widths is 
seen to drop significantly at ~ — 8A. Using this value as our 
EW{[0 n]} division gives a false positive to real detection 
ratio of 43.4 per cent. Decreasing to -5A gives a ratio of 52.7 
per cent, while increasing to -20A gives a ratio of 54.9 per 
cent. 

4.2 Spectral Properties of 2SLAQ LRGs 

Using the selection criteria outlined in the preceding sec- 
tion, we can determine the fraction of LRGs in each of the 
four categories. Table|5|shows the number of galaxies in each 
category, and what fraction of the total population they con- 
stitute. 

Errors on the stated fractions are calculated assuming 
Poisson statistics. 

The corrected fractions are calculated by subtracting 
the estimated false detection rate calculated in Section [4.11 
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Table 1. Passbands and side continuum bands used to measure the mock H<5 and mock [O n] equivalent widths. 



Index 


Line passband (A) 


Blue continuum sideband (A) 


Red continuum sideband(A) 


H<5 


4083.5 - 4122.25 


4041.6 - 4079.75 


4128.5 - 4161.0 


[O n] 


3713 - 3741 


3653 - 3713 


3741 - 3801 


H 7 


4319.75 - 4363.5 


4283.5 - 4319.75 


4367.25 - 4419.75 


mock H<5 


4230.5 - 4242.25 


4165.6 - 4199.75 


4245.5 - 4271.0 


mock [O II] 


3590.0 - 3617 


3480.0 - 3550.0 


3617.0 - 3641.0 



Table 2. Fraction of each spectral type and properties of their combined spectrum. 



Spectral Type 


Number 


Fraction (per cent) 


Corrected fraction (per cent) 


D4000 


EW{[0 ii]} (A) 


EW{H<5}(A) 


EW{H 7 } (A) 


Passive 


4612 


81 ±1 


88 ±1 


1.75 ±0.01 


-1.0±0.1 


-1.67 ±0.08 


-5.83 ±0.08 


k+a 


225 


3.9 ±0.3 


2.7± 0.2 


1.53 ±0.01 


-2.3 ±0.1 


3.7 ±0.1 


-1.7 ±0.01 


em+a 


101 


1.8 ±0.2 


1.2 ±0.1 


1.50 ±0.01 


-12.1 ±0.1 


3.9 ± 0.1 


-1.6±0.1 


cm 


759 


13.3 ±0.5 


8.6 ±0.4 


1.71 ±0.01 


-11.7 ±0.1 


-1.37 ±0.09 


-5.53 ± 0.08 




EW (A) 

Figure 2. Distribution of real H<5 (top) and [O Ii] (bottom) equiv- 
alent widths (solid line) and mock H<5 (top) and [O Ii] (bottom) 
equivalent widths (dashed line). The low signal to noise of the 
2SLAQ spectra is apparent from this figure as the mock equiva- 
lent width distributions cover a similar range to the real equiv- 
alent width distributions despite the fact that the there should 
be no strong features in the passband regions in the mock analy- 
sis. The dotted lines shows the equivalent width thresholds which 
define the classification scheme used throughout the paper. 



from the measured fraction. The false detection rates for the 
H<5 and [O n] line are transformed into rates for the k±a, 
em+a and em classifications by assuming that the false de- 
tection rate for a line is spread evenly over the classifications 
requiring it. For example, 30.1 per cent of H<5 detections in 



k+a and em+a LRGs are false detections, thus if they are 
spread evenly across both classifications the false detection 
rate in k+a and em+a LRGs (due to only H<5) is 20.8 and 
9.3 per cent respectively. In the case of the em+a classifica- 
tion we define the false detection rate to be sum of the false 
detection rate from both H5 and [O n]. 

The efficiency of the selection method in isolating pas- 
sively evolving ellipticals is again apparent from Table|5|with 
over 80 per cent of LRGs having the spectral features indica- 
tive of an old passively evolving stellar population. 

A number of example spectra for each of the spectral 
classes are shown in Figure [3] Spectral features discussed 
are marked. It can be seen that in the spectra with strong 
H<5, the other, higher order, Balmer lines are also observed. 

The combined spectra for each spectral type are pre- 
sented in Figure 0] The combined spectra are produced by 
first de-redshifting each spectrum, and rebinning to a com- 
mon dispersion. The spectra are first normalised using the 
mean flux in a 100A window (4000-4100A). The top and bot- 
tom 10 per cent of values at each pixel are then excluded. 
The mean of the remaining values at each pixel is then used 
to produce the combined spectrum. For each combined spec- 
trum the D4000 index and equivalent widths of the H5, H7, 
H/3 and [O 11] lines have been measured using the methods 
presented in Section ^. II with the results presented in Table 

m 

While it is generally accepted that H5 is an efficient in- 
dicator of a young stellar population, the presence of strong 
H<5 absorption should be accompanied by increased absorp- 
tion in the other Balmer lines. To test this we measure the 
equivalent width of the H7 line in the combined spectra, 
with the results presented in Table [5] H7 alone is utilised 
as the lower order Balmer lines (H/3, Ha) are not inside the 
wavelength coverage of many of the higher redshift spec- 
tra. It can be seen that for classifications which pass the 
EW{H5}> 2 A cut, k+a and em+a, the equivalent width of 
H7, is significantly greater in absorption than that found in 
the passive combined spectrum. It should be noted that the 
H7 equivalent widths quoted in Table|5|appear to be in emis- 
sion (i.e. negative) as emission in the nearby passbands used 
to estimate the continuum dominates moderately low levels 
of absorption in the line itself. The ratio of the Calcium H 
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3800 



4000 



4200 



4400 



Figure 3. Example spectra for the k+a (top three), em+a (mid- 
dle three) and em (bottom three) spectral classifications. Spectra 
have been converted to rest frame wavelengths, and smoothed by 
a Gaussian with FWHM=0.8A. Spectral features discussed in the 
text are marked. 



to K lines is also noticably greater in the k+a and em+a 
combined spectrum than in the Passive spectrum (Figure 
4), This can be attributed to increased absorption in the 
He line, which is coincident with the Calcium H line. Thus 
we can be satisfied that our selection based on H<5 alone 
has been effective in isolating LRGs with overall increased 
Balmer line absorption, and hence potentially younger stel- 
lar populations. 

Figure shows the D4000 distributions for the k+a 
(bottom), em+a (middle) and em (top) spectral classifica- 
tions. The whole 2SLAQ sample is designated by the solid 
line, while the D4000 distribution of only the passive LRGs is 
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Figure 4. Combined spectrum for the four spectral classifica- 
tions. Quantities measured on the spectra are presented in Table 
|3] Note the increase in the absorption strength of the higher order 
Balmer lines (H7,He, etc) in the the k/cm+a combined spectra. 



shown by the dashed line. The D4000 index is well known to 
be sensitive to recent star formation, with lower values indi- 
cating the increasing presence of a young A star population. 
As expected the D4000 indices for LRGs in the 2SLAQ sam- 
ple suspected to have recent star formation, (k+a & em+a) 
are noticably less than that for the passively evolving LRGs. 
The em+a LRGs also show significantly lower D4000 indices 
than the k+a LRGs suggesting that the em+a galaxies have 
experienced more recent, or ongoing given the presence of 
the [O n] line, star formation. Interestingly the em LRGs 
have D4000 indices which are broadly consistent with the 
passive LRGs. While a KS test shows that all three distri- 
butions are statistically inconsistent with the passive popu- 
lation, with p-values of 2.4xlCT 6 , 2.6xlCT ir and 6.6xlCT 22 
for the em, em+a and k+a LRGs respectively, a significant 
fraction of the em LRG population are found to have D4000 
indices as large as those found in the passive population. 
This is in stark contrast to the k+a and em+a LRG pop- 
ulations which are very rarely found at high D4000 index 
(i.e. D4000> 2). This would suggest that while recent star 
formation is unambigously responsible for k+a and em+a 
LRGs, the origins of em LRG population are not so clear. 
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Figure 5. Distribution of D4000 indices for k+a (bottom), em+a 
(middle) and em (top) LRGs. In each case the distribution is 
marked with a solid line, while the distribution of passive LRGs 
is marked with a dashed line. The k+a and em+a LRGs show sys- 
temically lower D4000 indices than the passive LRGs. This, along 
with their increased H5, confirms the presence of a younger stel- 
lar population in these LRGs. Significant numbered of em LRGs 
show D4000 strengths that are as large as those found in the pas- 
sive population, suggesting that a large fraction of em LRGs do 
not possess significant populations of young stars. 

4.3 Redshift Dependence 

To test the redshift dependence of the classification frac- 
tions, the 2SLAQ LRGs are divided into 2 subsamples be- 
tween 0.45 < z < 0.55 and 0.55 < z < 0.65. In order to 
increase the redshift range of the analysis a low redshift 
sample for comparison is provided by SDSS MAIN survey 
galaxies with 0.1 < z < 0.2. While the SDSS also conducted 
an LRG survey the selection criteria for SDSS LRGs is much 
stricter than for 2SLAQ, so strict that only ~700 2SLAQ 
LRGs would be selected by in the SDSS survey at z = 0.2. 
It is for this reason that we utilise galaxies from the SDSS 
MAIN survey for comparison. 

As the 2SLAQ LRGs were selected by an observed- 
frame colour selection some effort must be made to ensure 
our comparative sample from SDSS is consistent with the 
colour selections of 2SLAQ. Indeed effort must also be made 
to ensure that the colour selections are consistent within the 
redshift range covered by 2SLAQ itself. The simplest way to 
accomplish this is to use an independent rest frame selection 
in colour-magnitude space on both the 2SLAQ and SDSS 
MAIN galaxies. To determine the appropriate boundaries of 
this selection we present the M; and rest frame g—i evolution 



redshift 

Figure 6. Rest-frame g — i and Mi redshift evolution for the 
2SLAQ LRGs. Rest-frame quantities are produced via k+e cor- 
rections based on a Bruzual and Chariot (2003) model track. The 
model assume a single epoch of star formation for the LRG stel- 
lar population at high redshift, followed by passive evolution. The 
dashed line shows the cuts utilised in the text to produce a sam- 
ple of LRGs which are unbiased in colour and luminosity with 
redshift. 



of the 2SLAQ galaxies with redshift in Figure |5] Absolute 
magnitudes for the 2SLAQ galaxies are produced via the 
use of model k+e corrections based on Bruzual and Chariot 
(2003;henceforth BC03) models. The models assume sim- 
ple passive evolution of the stellar population from a single 
epoch of high redshift star formation. This sort of model has 
been found to approximate the colour evolution of LRGs to 
~ 0.1 mag (Wake et al. 2006), which is satisfactory for our 
purposes here. While ideally k-corrections could be derived 
directly from the spectra, the 2SLAQ spectra are not flux- 
calibrated, nor do they have the significant wavelength cov- 
erage required for such calculations. Significant redshift evo- 
lution is shown in both the g — i colour and Mi of the LRGs 
in Figure [U The Mi evolution is a result of the iDev < 19-8 
magnitude limit of the 2SLAQ survey, while the increase in 
the distribution of g — i colours to bluer values is a result of 
the colour selection becoming 'looser' at higher redshifts (as 
can be seen to some extent in Figure It can be seen that 
a cut at g — i > 0.8 (as shown by the dashed line) will pro- 
duce a sample which is reasonably consistent in colour with 
redshift. Similarly a cut at Mi < —22 will produce a sam- 
ple which is close to volume-limited. It is on this cut-down 
'homogeneous' sample of 2SLAQ LRGs that the following 
analysis is performed. 
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In a similar fashion we select galaxies from the SDSS 
MAIN survey which have rest-frame g — i > 0.9 and Mi < 
—22.0. As the SDSS spectra have a wide wavelength cover- 
age, and good spectrophotometric calibration k corrections 
can be derived directly from the spectra. This has been done 
for the DR4 MAIN galaxies by Blanton et al. (2005) and can 
be found in the NYU SDSS value-added catalogue (VAC) 
available at http://wassup.physics.nyu.edu/vagc/ K 
corrections in the NYU VAC are produced via the use of the 
kcorrect software (v3.4) described in Blanton et al. (2003). 
As the VAC only provides k-corrected photometry the Mj 
cut is adjusted slighty to -22.2 to take into account of the 
Mi evolution from z = to z — 0.15 predicted by the BC03 
models used to k+e correct the 2SLAQ LRGs. 

The results of this analysis are presented in Figure |7| 
and Table |H] In each subsample from 2SLAQ, and for the 
low-z SDSS galaxies, we calculate the number of each spec- 
tral classification present, and the fraction of the bin that 
this number constitutes. Errors are again calculated assum- 
ing Poisson statistics. Since the LRGs are not uniformly dis- 
tributed inside each bin the bin centres are taken to be the 
median redshift of the all the LRGs in a given bin. 

In order to correct for the occurrence of false detections 
we again apply corrections to the fractions estimated using 
the results of the 'mock' line analysis from Section [4.11 By 
dividing the 'mock' equivalent width measures on H<5 and 
[O il] into the same redshift bins as the real analysis we 
produce corrections which are indicative of the underlying 
signal to noise distribution in each bin. This has the ad- 
vantage of also eliminating any possible bias that may arise 
from varying signal to noise across the redshift range. 

The raw and corrected fractions are presented in Figure 
|7|with the values shown in Table [3] 

The effect of the g - i > 0.9 and Mi < -22 cuts on the 
number of k+a and em+a LRGs is quite dramatic. Only 79 
out of the original 225 k+a's and even more dramatically 
only 23 out of the original 101 em+a LRGs pass the cuts. 
However this is not all that surprising given the strong corre- 
spondence between k+a and em+a features and low D4000, 
which can be taken as an effective proxy for g — i colour. 

It can be seen from Figure |7| that only the k+a LRG 
fraction shows a clear trend of increasing with redshift. The 
em LRG fraction shows a strong increase between the SDSS 
z = 0.17 sample and 2SLAQ LRGs at z = 0.49, but then 
shows a significant drop in the highest redshift bin. Little 
can be said about the em+a LRG fraction as the numbers 
in the 2SLAQ bins are too low. 

In order to quantify the redshift dependence of the k+a 
fraction in Figure[7|a (l + z) n fit has been performed on the 
corrected fractions. Fitting is performed via weighted linear 
regression in log(z) vs log(fraction) space. The fit is shown 
by the dashed line in Figure[7| The k+a fraction is found to 
obey a (1 + z ) 2 - 8±0 - 7 relation. 



5 DISCUSSION 

5.1 Comparison with previous work 

In recent years many authors have looked for recent and/or 
ongoing star formation amongst samples of early-type galax- 
ies. Fukugita et al. (2004) performed a spectroscopic study 
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Figure 7. Redshift dependence for the raw (crosses) and cor- 
rected (squares) fraction of each classification with redshift. For 
comparison the fraction of each classification in a sample of low 
redshift (z ~ 0.15) SDSS MAIN galaxies is shown. The SDSS 
MAIN galaxies are required to have the same luminosity and 
colour range as the 2SLAQ LRGs by placing simple cuts in rest- 
frame g — i and Mi. While the em and em+a fractions do not 
show any overall trend with redshift, the k+a fraction shows a 
strong increase. Fitting a (1 + z) n relation to the k+a fraction 
increases gives n = 2.8 + 0.7. The fit is shown by the dashed line. 

of 420 E/S0 z ~ 0.12 galaxies selected from SDSS. They 
found evidence of ongoing star formation, in the form of Ha 
emission, in 19 (4.5 per cent) and k+a like spectral prop- 
erties in 15 (3.6 per cent). Yi et al. (2005) also identified 
early-types with recent star formation in SDSS, this time us- 
ing GALEX Near-UV and SDSS r-band photometry. ~ 15 
per cent of their sample of 39 bright (M r < —22), early-type, 
z < 0.13 galaxies was found to have evidence of recent star 
formation via the use of NUV-r colour. Moving to higher 
redshifts Doherty et al. (2005) and Le Borgne et al. (2005) 
both find significant fractions (~ 30, and 50 per cent, re- 
spectively) of k+a's amongst massive galaxies at z ~ 1. 

A significant fraction of 2SLAQ LRGs fall into the em 
category (6.8 ± 0.3 per cent), although this result is not 
completely unexpected as a study by Willis et al. (2002) 
found over 20 per cent of their sample of luminous field 
early type galaxies had significant [O n] emission. Similarly 
Eisenstein et al. (2003) found that ~ 10 per cent of their 
sample of massive bulge dominated galaxies had emission 
line components. 

However it is difficult to make concrete comparisons be- 
tween these results and those presented here. This is because 
the fraction of early-types with young stellar populations is 
very sensitive to the sample selection. This can even be seen 
here, with significant differences in the fraction of k+a and 
em+a LRGs between the whole sample in Section 14.21 and 
cut-down sample used in Section 14.31 Most previous stud- 
ies, including those discussed above, use a selection based 
on morphological information. This is distinctly different 
from the colour based selection used in 2SLAQ survey. It 
should be expected that morphology based studies of early- 
types will find greater fractions of galaxies with young stel- 
lar populations than we find here as the colour selection is 
strongly biased against the selection of galaxies with large 
populations of young stars. Thus while the absolute frac- 
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Table 3. Redshift dependence of classification fractions. Quoted errors are calculated assuming Poisson statistics. The fractions are 
calculated by first producing a truncated sample which is uniform in colour and luminosity with redshift, and then correcting the 
fractions using the same analysis as described in Section 14.11 The low redshift z ~ 0.17 point is measured using SDSS MAIN galaxies 
with the same rest-frame colours and luminosities as the 2SLAQ LRGs. 
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tions of k±a and em±a LRGs quoted here are not directly 
comparable to other work on massive early-types, the ob- 
served increase with redshift should be consistent between 
the studies. 

This notion is not only supported qualitatively by the 
jump in the number of 'young' early-types between z ~ 0.1 
and z ~ 1 but also quantitatively by a comparison of the 
measured redshift dependence of the k±a fraction here and 
in the study of Le Borgne et al. (2005). Le Borgne et al. 
(2005) performed a study of the spectral properties of mas- 
sive (M> 10 10 02 Mq) galaxies selected from the Gemini 
Deep Deep Survey (GDDS; see Abraham et al. 2004) and 
SDSS MAIN galaxy sample. While they did not base their 
selection on morphology or colour, and hence found a sig- 
nificantly greater number of k±a than we find, they found 
an evolution of (l±z) 2 ' 5±0,7 across a redshift range of the 
0.1 < 2 < 1.2. Le Borgne et al. (2005) note that this value 
is very close to the predicted evolution of the galaxy merger 
rate in ACDM semi-analytic models of (1 + z) 3 ' 2 (Le Fevre 
et al. 2000). Interestingly the fraction of k+a and em±a 
LRGs found here is very close to 4 ± 1 per cent fraction of 
massive galaxies in close pairs found by Bell et al. (2006). 
Thus if k±a and/or em±a activity is triggered by merg- 
ers, the fractions and evolutionary trends we find here are 
reasonably consistent with other merger indicators in the 
literature. 



5.2 H<5 strong LRGs 

If we accept the interpretation that the observed H(5 absorp- 
tion in the k±a and em±a LRGs is a result of recent star 
formation, which seems likely given the weaker D4000 in- 
dices and observed absorption in the higher order Balmer 
lines (see Section 14.20 . then a small, but significant, num- 
ber of the LRG population (~ 1 per cent) has been forming 
stars within the last 2 Gyr. This is an important result as it 
shows that evolution is occurring in the massive early-type 
population at redshifts less than 0.7. 

To try and quantify the level of star formation required 
to produce the observed EW{H<5} we turn to the spectral 
synthesis models of Bruzual & Chariot (2003). A distinct 
advantage of using these models is that they produce output 
measurements of the indices discussed here (D4000 and H<5) 
measured in same fashion as outlined in Section|3] All of the 
models presented in this Section assume a Salpeter (1955) 
IMF, solar metallicity, and no dust reddening. 

While constant star formation in a galaxy will produce 
the required level of H<5, this scenario has two major prob- 
lems. Firstly k±a LRGs show no [O n] emission, suggesting 
that there is no ongoing star formation occurring, however 



it is possible that the [O n] could be missing due to pref- 
erential dust obscuration or even aperture effects. Secondly, 
models show that galaxies with continuous star formation 
will mostly be outside of the LRG colour selection. Figure|H] 
shows model tracks for two constant star formation models 
and a purely passive evolution model. In each of the con- 
stant star formation models a large fraction of the stellar 
population is formed at high redshift [z ~ 3), with low lev- 
els of constant star formation responsible for 10 and 5 per 
cent of the stars at age 10 Gyr. The 10 per cent model track 
falls well outside of the 2SLAQ colour selection, while the 
5 per cent model track just enters the 2SLAQ selection at 
z > 0.55. The typical EW{H5} of the model spectra in the 
redshift range probed by 2SLAQ is ~ 2.9A and ~ 1.6A for 
the 10 per cent and 5 per cent model respectively. Thus 
the 10 per cent model spectra would be selected in our k+a 
sample, but not in 2SLAQ and the 5 per cent model would 
be selected in 2SLAQ, but not as a k±a. As the 5 per cent 
model only just falls into the 2SLAQ selection we consider 
it unlikely that continuous star forming models can explain 
the k±a LRGs in 2SLAQ. 

Thus we confine ourselves to models in which a sig- 
nificant episode of star formation has occurred and subse- 
quently ceased. In order to simplify the models we only con- 
sider scenarios in which an instantaneous burst of star for- 
mation is superimposed on an old stellar population. While 
it is likely that star formation histories with longer periods 
of star formation are responsible for the k/em+a popula- 
tion, the tracks taken by such models after star formation 
has ceased are nearly indistinguishable from the instanta- 
neous burst models. The only quantitative difference is that 
longer star burst models require that a larger fraction of the 
LRGs mass be consumed in star formation to achieve the 
same level of H<5 absorption. This can be seen from Figure|U] 
which shows the relationship between burst length, and peak 
EW{H<5} for models which consume 2, 5, and 10 per cent of 
galactic stellar mass in new star formation. In all cases the 
models superimpose a burst of constant star formation on 
an old (7.4 Gyr) stellar population. The peak EW{H<5} is 
found to not only be dependent on the size of the burst (i.e. 
the total number of new stars created) but is also strongly 
dependent on the period over which the burst occurs. Indeed 
for bursts in the range 500-2000 Myr the peak EW{HW} is 
dominated by the star formation rate, not the total size of 
the burst. This can be seen if we take an arbitrary peak 
EW{H<5} and look at the star formation rate (in per cent of 
galactic stellar mass per year) required for the three burst 
sizes presented in Figure E] Using EW{H<5}=4A as an ex- 
ample, we estimate star formation rates for the 2 per cent, 
5 per cent, and 10 per cent models of 6.6 x 10~ 9 , 4.2 x 10~ 9 
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Figure 8. Model g — r, r — i tracks for two constant star for- 
mation models (crosses and asterisks), and a passive evolution 
model (triangles). The constant star formation models assume 
most of their stars are formed at high redshift, with a low level 
of constant star formation occurring since. The crosses represent 
a model in which 90 percent of the stellar population at age 10 
Gyr was formed at high redshift, the asterisks represent a model 
where 95 per cent were formed at high redshift. The dashed line 
shows the Sample 9 2SLAQ selection described in Sectional The 
passive evolution model assumes that all of the stars are formed at 
high redshift, and have experienced only passive evolution since. 
The circles indicate redshift of z = 0.4 and the tracks mark the 
evolution from z=0.8 (near the upper left border) to z = 0.4 at 
Az = 0.025 intervals. Assuming a formation age of z = 3 for 
the model, the typical EW{H<5} in the redshift range probed by 
2SLAQ is ~ 2.9A and ~ 1.6A for the 10 per cent and 5 per cent 
model respectively. Thus it can be seen that any galaxy experi- 
encing a level of constant star formation high enough to pass our 
EW{Hi5}> 2A cut would not be selected by the 2SLAQ colour 
selection. 
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Figure 9. Effect of changing the length of the star burst on the 
peak EW{H<5}. Bursts consuming 2 (crosses), 5 (circles), and 10 
(asterisks) per cent of galactic stellar mass are shown. The peak 
EW{H<5} is seen to decrease with increasing the burst length. It is 
likely that the star formation rate, not the total size of the burst, 
is the dominating factor in determining the peak EW{H<5}. 

plausible amounts, can only move the tracks shown in Fig- 
ure GUI by at most 0.1 in D4000 (Balogh et al. 1999, Pracy 
et al. 2005), not nearly enough to explain the reddest H<5 
strong LRGs. Curiously the fraction of red H<5 strong LRGs 
in 2SLAQ is exactly equal to the predicted number of false 
positives in the sample. Given this it would seem probable 
that most, if not all, of the red H<5 strong LRGs in this study 
are simply false detections. If this is true than the origin of 
the k+a and em+a LRG population can be most readily ex- 
plained by recent bursts of star formation, on the scale of 1 
per cent of galactic stellar mass, in otherwise passive LRGs. 



and 5.0 x 10 9 per cent per year, respectively. It can be seen 
that these values are in relatively good agreement. 

Figure ESI plots the track taken in EW{H5}- D4000 
space by models in which star bursts consuming 1,2, and 5 
per cent of stellar mass have occurred. In each case the model 
assumes an instantaneous burst occurs in a passively evolv- 
ing stellar population which formed 7.4 Gyr before burst. 
Also shown in Figure 1101 is the number density of LRGs 
plotted in grey scale. As expected a slight trend of increas- 
ing EW{H5} with decreasing D4000 is observed. Reasonable 
agreement is found between the H<5 strong LRGs and the 
model tracks. Given that we have only presented tracks for 
a limited range of models it should be possible to produce 
models in a similar vein to those presented which can match 
any observed EW{H5}. However this is not the case for all 
values of D4000. The dotted vertical line at D4000=1.65 is 
used to crudely divide between those H5 strong LRGs which 
roughly agree with the models and those which don't. Of 
the 326 H<5 strong LRGs in 2SLAQ 98 (30.1 per cent) have 
D4000> 1.65. Many of the previous E+A/H5-strong studies 
have also identified a population of red H<5 strong galaxies 
with photometric properties that are too red to be consis- 
tent with models (Couch & Sharpies 1987; Blake et al. 2004; 
Poggianti et al. 1999; Balogh et al. 2005). While dust and 
super solar metallicities can be called upon to redden the 
model tracks, and hence explain the discrepancies observed, 
it has been repeatedly shown that these effects, when used in 



5.3 em LRGs 

The origin of the em LRGs is more ambiguous as the ob- 
served [O n] emission could originate from either ongoing 
star formation or AGN. Typically the nature of emission 
lines in optical spectra can be determined via the use of line 
ratios ( e.g. Baldwin et al. 1981), or even more unambigously 
by their radio properties. Unfortunately neither of these di- 
agnostics are readily applicable here. The traditional opti- 
cal line ratios used for classifying AGN/SF are H/3, [OIII] 
5007A, and Ha, all of which are have wavelengths too red 
to appear inside the wavelength coverage of 2SLAQ spectra. 
Sadler et al. (2006) have performed a crossmatching of the 
2SLAQ LRGs to radio sources in both the FIRST and NVSS 
catalogues. Of the 13784 2SLAQ LRGs with confident red- 
shifts only 378 have corresponding radio sources in FIRST, 
this is despite FIRST covering 96.5 per cent of the 2SLAQ 
survey area. Of the 759 em LRGs found in the sample pre- 
sented here, only 21 (2.8 per cent) have a FIRST radio detec- 
tion within 5" , very close to the overall radio detection rate 
of 2.7 per cent. Similarly Sadler et al. (2006) find that the 
fraction of radio detections with [O n] emission is the same 
as the fraction in the whole 2SLAQ sample. However while 
strong radio continuum emission can unambigously identify 
an AGN, the converse, i.e. lack of radio emission constitutes 
a lack of AGN activity, is not true. Thus while we can be 
confident that at least 2.8 per cent of em LRGs have AGN, 
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Figure 10. Number density of LRGs in the D4000-EW{H<5} 
plane. The plane is broken down into 0.1x1 cells, the grey scale 
representing the logarithm of the number of LRGs in the cell. 
The scale goes from light to dark, with the darkest spot equal to 
316 (10 2,5 ) LRGs. Also plotted are model tracks for a series of 
models assuming bursts of star formations consuming 5 (top), 2, 
and 1 (bottom) per cent of stellar mass in new star formation. 
Typical errors on both indices are shown in the top right corner. 
A significant number of LRGs have EW{H5} and large D4000 
indices which cannot be explained by these simple models. These 
are analogous to red H<5-strong galaxies found in the literature 
(e.g. Couch & Sharpies 1987). 

little can be said about the remaining 97.2 per cent from 
available radio data. 

However despite the lack of data with which to accu- 
rately classify each em LRGs, we can muse about the origins 
of em LRGs. A strong piece of evidence against a star form- 
ing origin for all em LRGs is the lack of correlation between 
[O n] emission and lower D4000 indices. The D4000 index 
is very sensitive to any population of young stars, as seen in 
the case of the k+a and em+a LRGs in Figure^] Thus un- 
less we are observing an LRG during the very short period 
of time at the very beginning of a star formation episode be- 
fore substantial numbers of young stars have been formed, 
their should be a clear relation between the D4000 index 
and current star formation. However if the [O n] emission 
originates from AGN activity this should have no discern- 
able effect on the D4000 index. Given this it is unlikely that 
the large number of em LRGs with relatively large D4000 
indices have star-forming origins. However given the large 
spread in the D4000 index seen in the passive LRG popula- 
tion (Figure |5J it is impossible to say whether the em LRGs 
with low D4000 index are actively star forming. 

5.4 Recent Star Formation in LRGs and the 
ACDM Paradigm 

Assuming that recent and or ongoing star formation is the 
cause of the observed spectral characteristics of the k/em+a 



LRGs, which is likely given the arguments of the previous 
sections, a number of well-known phenomena can be called 
upon as potential candidates for origin of star formation in 
LRGs. 

Given that LRGs are akin to massive early-type galaxies 
a sizable fraction will reside in groups or clusters, making a 
small, but significant, number of them bright cluster or even 
cD galaxies. Thus it is plausible that these central cluster 
LRGs may be at the focus of a cooling flow. The presence 
of optical emission lines in 20-30% of all brightest cluster 
galaxies in X-ray selected samples has been confirmed in 
a number of studies (Donahue et al. 1992; Crawford et al. 
1999) and an extreme population of galaxies with massive 
starbursts (50-200 M©yr _1 ) has been identified (Allen 1995; 
Crawford et al. 1999). These galaxies could be the progen- 
itors of the k/em+a LRGs if the starburst is truncated al- 
though they would be preferentially found in the brightest 
LRGs. 

Another possibility is that the bursts of star formation 
are being created as a result of LRG-galaxy mergers or inter- 
actions. This ties in well with the idea that LRG evolution is 
dominated by hierarchical merging. Indeed in this model one 
would naively expect (and indeed it has been shown from 
N body simulations) that the frequency of mergers would 
decrease significantly with increasing time (decreasing red- 
shift), which is qualitatively in line with the results demon- 
strated here. 

However it is difficult to make quantitative assertions 
about the origin of k/em+a LRGs as it is impossible to de- 
termine which of these hypotheses is responsible, if at all. 
Another complication is the potential for differences between 
the k+a and em+a classes. Previous studies of traditional 
k+a's and em+a analogues at low redshift have found sys- 
tematic differences between the two classes, suggesting dif- 
ferent formation mechanisms (Balogh et al. 2005). 

Regardless of the physical nature of the star formation 
in LRGs, it is reasonably safe to assume that the star form- 
ing gas is coming from an external origin. This is in agree- 
ment with the core principles of ACDM cosmology. Results 
from the current generation of hydrodynamic N-body cos- 
mological simulations suggest that the most massive galax- 
ies increase in mass significantly from z < 1 (De Lucia et al. 
2005), with the bulk of this increase coming from accretion of 
intra-cluster gas, or via mergers with smaller galaxies (Mu- 
rali et al. 2002) . By consideration of the expected timescale 
of the k+a signature we can crudely estimate the number 
of LRGs affected by k+a activity since z = 0.8. Converting 
the redshifts of the bins used in Section 14. . SI into look-back 
time results in a 0.36 x ^°- 85±0 - 2 dependence, where t is the 
look-back time in Gyr. By simply integrating we can use this 
dependence to estimate the fraction of LRGs which undergo 
a k+a phase between < z < 0.8 (or < t < 7) as 8 per 
cent. For comparative purposes we consider both the mod- 
els of De Lucia et al. (2005) and the observational evidence 
presented by Bell et al. (2006). De Lucia et al. (2005) pre- 
dict that ~ 50 per cent of massive early-types gain 50 per 
cent of their stellar mass between < z < 0.8. Similarly Bell 
et al. (2006) predict from there measure of the close pair 
fraction that ~ 20 per cent of massive galaxies will undergo 
a major merger between < z < 0.8. Thus even if k+a, or 
em+a, activity was associated with a major merger event 
that doubled the mass of the progenitor LRG, which is un- 
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likely given the models presented in Section 5.2, this would 
only account for ~ 16 per cent of the required mass growth 
in massive early-types in the models, and only 40 per cent 
of the predicted number of mergers in massive galaxies from 
observations. It is clear from this crude comparison that the 
level of k+a and em+a activity is not sufficient to explain 
the bulk of the mass growth in massive early types. This is 
in line with recent results which suggest that red or 'dry' 
mergers are responsible for the growth of massive early-type 
galaxies (van Dokkum et al. 2005; Bell et al. 2005). 



6 CONCLUSION 

We have determined the recent star formation histories of a 
sample of 5,697 LRGs based on the equivalent widths of the 
H<5 and [O n] lines. While the majority (>80 per cent) show 
the spectral properties of an old, passively evolving, stellar 
population, a significant number of LRGs show evidence for 
recent and/or ongoing star formation in the form of k+a 
(2.7 per cent), em+a (1.2 per cent) or em LRGs (8.6 per 
cent). By dividing the sample into 2 redshift subsamples 
from 0.45 < z < 0.55 and 0.55 < z < 0.65, and comparing 
to a z ~ 0.15 sample selected from SDSS, it is observed 
that the fraction of k+a LRGs increases with redshift as 
(l + z) 2 - 8±0 - 7 . 

Spectral synthesis models, utilising the code of Bruzual 
& Chariot (2003), suggest that the k/em+a LRGs could 
originate from passive LRGs which undergo a starburst in 
which material equivalent to ~ 1 per cent of the stellar mass 
is consumed in new star formation. 

Several origins for k+a and em+a LRGs are considered, 
including cooling flows and mergers, however identification 
of the exact formation mechanism, or even if k+a and em+a 
LRGs share a common origin, is not possible with the cur- 
rent data. It is clear that k+a and em+a activity represents 
a signpost of recent evolution in LRGs, however by consid- 
ering the life time of the k+a/em+a phase we estimate that 
only 8 per cent of LRGs will experience a k+a or em+a 
phase between z = 0.8 and the present. By comparing this 
with the semi-analytic galaxy formation models of De Lucia 
et al. (2005) we conclude that k+a/em+a activity cannot 
be responsible for the predicted mass growth of massive el- 
lipticals since z = 0.8. 
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